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One of the most pressing issues in cosmology is whether general relativity (GR) plus a dark 
sector is the underlying physical theory or whether a modified gravity model is needed. Upcoming 
dark energy experiments designed to probe dark energy with multiple methods can address this 
question by comparing the results of the different methods in constraining dark energy parameters. 
Disagreement would signal the breakdown of the assumed model (GR plus dark energy). We study 
the power of this consistency test by projecting constraints in the wo—Wa plane from the four difi'erent 
techniques of the Dark Energy Survey in the event that the underlying true model is modified gravity. 
We find that the standard technique of looking for overlap has some shortcomings, and we propose 
an alternative, more powerful Multi- dimensional Consistency Test. We introduce the methodology 
for projecting whether a given experiment will be able to use this test to distinguish a modified 
gravity model from GR. 

PACS numbers: 95.35.+d; 95.85.Pw 



I. INTRODUCTION 



General relativity (GR) is currently a bad fit to cosmological data unless a new substance, so-called dark energy, 
is invoked. If GR really is an incomplete or incorrect theory and we are tasked with identifying the correct model, 
a major hurdle will be determining how to confront upcoming data sets in the absence of a well-understood model. 
What new parameters should be introduced and fit for when, e.g., data on weak gravitational lensing or galaxy clusters 
are analyzed? Several authors have addressed this question [iHa], and it has recently become possible to test GR 
using survey data 

Here we address a slightly less ambitious question: using multiple cosmological probes, how can we determine 
whether cosmic acceleration is driven by dark energy or modified gravity (MG)? One approach is to analyze the data 
assuming that GR is correct and see whether the constraints on dark energy parameters from different probes overlap 
[lol [Tl|. Non-overlapping constraints would be a strong signal that the underlying parameterization is wrong; i.e, 
that GR-|-dark energy cannot account for the data and that a modified theory of gravity is called for. A similar 
approach is to look at parameter constraints coming from separate dynamical effects such as the cosmic expansion or 
perturbation growth jl2l | . Here we explore the former method in depth in the context of a concrete example. 

Ishak et al. showed that, in pri ncip le, non-overlapping dark energy parameter constraints obtained from multiple 
experiments is a signature of MG In particular, they found that dark energy parameters obtained from a space- 
based supernova survey and a space-based weak lensing survey will not agree if the Universe is in fact described by the 
Dvali-Gabadadze-Porrati (DGP) braneworld model We reexamine this general method with our own example, 
assuming that the universe is governed by a toy MG model and considering projections from the upcoming Dark 
Energy Survey (DES). We present the projected constraints from all four DES probes in the plane of dark energy 
parameters wq and Wa, where the dark energy equation of state is assumed to be w = wq + ^0(1 — a) and a is the scale 
factor of the universe. This straightforward plot is not the most powerful way to combine probes, so we introduce a 
more quantitative formalism that should be useful for future attempts in this direction. The formalism assigns a 
for the combined probes which can be interpreted in the usual fashion so that a "bad" corresponds to disagreement 
among the probes, and therefore a quantitative assessment of how well the model of GR-|-dark energy works. 

Section II discusses modified gravity models in general and details the modified gravity model we adopt as our 
working example. Section HI then presents the DES projections in the {'Wo,'Wa) plane along with a description of 
the shortcomings of this approach. In Section IV, we present a more quantitative approach (see also [l3|). which we 
call the Multi- dimensional Consistency Test (MCT), illustrate how to obtain MCT projections, deal with the issue of 
degenerate directions, and finally conclude by applying this formalism to DES for the model under study. 
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II. PERTURBATIONS IN MODIFIED GRAVITY 

The metric in the class of modified gravity models we consider retains its standard GR form 

ds^ = - {1 + 2^) dt^ + {1 + 2<i>) dx^ (1) 

where a is the cosmic scale factor and $ and ^I* are the scalar gravitational potentials. Hu & Sawicki Q proposed 
introducing two functions which parameterize deviations from GR: 

$ + * 



where is the background matter density. In GR, the two scalar potentials are equal and opposite and Poisson's 
equations holds, so both / and g vanish. In modified gravity models, though, both / and g can be non-vanishing and 
are functions of both wavenumber k and scale factor a. 

Other parameterizations can be rewritten in terms of / and g. For example, instead of g, ij = —$/'!' is often used, 
so that g = (?7 — l)/(77 + 1). Zhao ct al. Q introduce a function which governs the peculiar velocity potential 

(3) 



inGprnO^S 

In terms of / and g, 

so ^ = 1 in standard GR. Another common way of parameterizing the effects of modified gravity was introduced by 
Linder P, Q and treats the growth factor of matter perturbations 

'^-n^iar (5) 
din a 

where rim(a) = rtm.o/[H{a)/HQ]'^. For a wide variety of models, taking 7 to be a constant works well and is an 
appealing way to confront data. In the context of the Hu & Sawicki formalism, 7 is not another new parameter, but 
rather is governed by the growth equation which becomes 



7-1 



dlnH 

1-27 — 37 + 2 

dm a 



, o ^ 3 3(1 - g) 



where H{a) is the expansion rate. In standard GR, one finds that 7 ~ 0.55. A choice of 7 corresponds to a choice of 
/ and g. Therefore, one cannot choose all three independently. 

In this work, we choose as the underlying "true" model a toy model in which 7 = 0.68 and / = 0. This is 
approximately the case in the Dvali-Gabadadze-Porrati (DGP) braneworld model [l|, lij, though we will not be using 
that model specifically and only mention it here as a motivation. The background expansion in DGP is different from 
that of ACDM, but we have chosen our model to have w = —1. Thus, the only observable differences between our toy 
model and ACDM will enter via the growth function, which determines the normalization of the linear matter power 
spectrum as a function of redshift. This is a mild modification to GR since scale-dependent growth modifications are 
indeed possible, and it has been shown that scale-dependent MG will be easier to detect 15|. In order to maintain 
consistency, our choices of 7 and / determine g via Fig. [1] shows the resulting g{z); note that constant 7 is not 
consistent with a constant g model. It is worth noting here that the MG model will produce more structure at early 
times than ACDM for fixed fluctuation amplitude today as- This is because structure grows more slowly in the MG 
model, therefore the Universe was clumpier (closer to its present state) in the past. 

III. DES PROJECTIONS FOR THE DARK ENERGY EQUATION OF STATE 

A. Parameter Misestimation 



For concreteness, we focus on the upcoming Dark Energy Survey, which will probe dark energy using thousands 
of Type la supernovae (SN), many high redshift clusters (CL), the Baryon Acoustic Oscillation (BAO) scale, and 
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FIG. 1: The post-Friedmann parameter (; as a function of redshift. Standard general relativity corresponds to g = or 7 ~ 0.55 
(dashed curve). The fiducial value used in this paper is shown by the solid curve corresponding to a model in which 7 = 0.68. 



the cosmic shear signal extracted from weakly lensed shapes of millions of background galaxies (WL). There is a 
wcU-cstablishcd formalism, the Fisher matrix approach, for projecting constraints from experiments such as DES. 
Ordinarily, one forms the Fisher matrix from which contours representing the 68% confidence region in the (wq^Wo) 
plane, say, can be drawn. These contours are drawn centered on the assumed underlying model. The Fisher matrix 
formalism is valid when the joint likelihood function of the cosmological parameters is a gaussian. 

Here we are after something a little different. We want to determine not only how large the error contours will be, 
but also where they would be centered in the event that an incorrect model is used to analyze the data. Our assumed 
true model is a toy modified gravity model with 7 = 0.68 and the same cosmic expansion as ACDM. We want to 
know what answer an analyst would get if s/he fit for dark energy parameters while assuming that GR was correct 
(7 = 0.55). What values of uig and Wa would be obtained? There is a simple extension [15[ of the Fisher formalism 
which provides an answer to this question. It involves three steps: 

1. Calculate the difference in the quantity to be measured (e.g., a power spectrum) in the true model and in the 
fitted model. Call this AP^ where subscript i labels the bins in which it is measured 

2. Calculate the Fisher matrix for the parameters Aq to be fit to the data. For a single experiment. 



F., = E(Cov ), — — 



(7) 



where Pi is the observed quantity in bin i, and (Cov)y is the covariance matrix for bins i and j. The covariance 
should be calculated using the cosmological model assumed to be true, while the derivatives should be calculated 
using the model we will fit. Priors may be added to the Fisher matrix if they come from probes which are robust 
in the case of MG (e.g. a prior on the Hubble constant, K). 



3. To first order in APj, the parameter Aq will be mis-estimated, or biased, by an amount 



AA„ 



/3 i3 



)..^AP. 



(8) 



where F^fj includes any priors. 
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FIG. 2: LEFT: Forecasted la constraints on dark energy parameters from the DES probes, including only statistical errors 
and assuming ACDM as the true model. From the largest to the smallest ellipse, the probes considered are baryon acoustic 
oscillations (black), supernovae (green), cluster counts (magenta), and weak lensing (blue). Each constraint is combined with 
a prior expected from Planck CMB measurements; additionally, the supernovae constraint includes an 8% prior on Ho- 

RIGHT: Same as LEFT but now the true model is assumed to be our toy modified gravity model with 7 = 0.68. 
Shown are the forecasted constraints when we incorrectly attempt to fit a GR+dark energy model to the data. The center of 
the weak lensing ellipse has moved to {wo, «;a) = (-l.l, 0.47) while the cluster counts ellipse has moved to (-1.19, 0.90). The 
probes are seemingly consistent, but we discuss the problems with this interpretation in Section [III Bl 



Our first task then is to determine the expected values of the measurements for the four probes in the assumed modi- 
fied gravity model and compare those to the predictions in standard GR+Dark energy. We consider a set of 8 standard 
cosmological parameters with fiducial values {wq, Wa, ^de, ^k,h, O;,, ns,as} = { — 1, 0, 0.73, 0, 0.72, 0.046, 1, 0.8} where 
Ofc is the curvature density, h is the Hubble constant in units of 100 km/s/Mpc, fif, is the baryon density, ris is the 
slope of the primordial spectrum, and erg normalizes the matter power spectrum at z = 0. For each probe, we then 
compute the constraints including projected priors from the Planck satellite [see e.g. [l^. We include only statistical 
errors in the projections for each experiment, therefore our parameter constraints will be optimistic but sufficient for 
our goal, which is to compare methods of testing GR. 

For two probes, supernovae and BAO, the answer is simple: these probes are sensitive only to background geometry 
which is assumed identical in our MG and GR models, so the predictions for the distance moduli (from supernovae) 
and correlation function peak (due to BAO) are identical to standard GR and AP = 0. The projected contours 
therefore are centered on the point in parameter space corresponding to the fiducial values. The only work that needs 
to be done is to determine the Fisher matrix which delineates the allowed region. This has been done before; here we 
simply reproduce these results, shown projected onto the {wo,Wa) plane in Figured] The CMB is mostly insensitive to 
our choice of MG since 7 only determines structure growth in the late Universe. The CMB power spectrum is in fact 
affected by gravity modifications via the late Integrated Sachs- Wolfe effect [13, ll^l ^-nd gravitational lensing, but we 
ignore these effects, which should only reduce our sensitivity to MG. Our Planck prior is therefore unchanged between 
the GR and MG cases. Only the weak lensing and cluster predictions are significantly changed when comparing GR 
to our toy MG model. Details on these calculations and Fisher matrix calculations for all probes are provided in the 
appendix. 
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B. Initial Results and Shortcomings 

The right side of Figure [2] suggests that if the underlying model were modified gravity with 7 = 0.68, then the 
allowed "la" regions from the 4 probes of DES are expected to overlap in the {wo,Wa) plane. They would show some 
disagreement but certainly not enough to reject the notion that dark energy and GR are the correct description of 
nature. This is a somewhat disappointing conclusion as the two models make very different predictions, so it behooves 
us to re-examine the proposed methodology. Perhaps the models can be distinguished with a more powerful statistical 
probe. 

The ambiguous conclusion stems from the approach of plotting the constraints from the four probes in the {wq, Wa) 
plane and seeing if they overlap. There are three shortcomings of this approach. First, the approach is not quantitative; 
it is visually appealing (do the contours overlap?) but does not yield a statistical conclusion. Second, each of the 4 
constraints are obtained with the Planck prior added in. This prior therefore is used multiple times. To the extent 
that one can obtain quantitative conclusions from observing the (non-)overlap of the contours, those conclusions will 
be incorrect due to the redundant information. Finally, the allowed regions are actually 8 dimensional. It is very 
possible that the allowed regions do not overlap in this 8D space but do in the projections onto the 2D subspace. 
Ishak et al. [ll| did find significant disagreement in the 2D subspace; however, compared to this work, their test-case 
considered more futuristic experiments (space-based supernovae and weak lensing) and a MG model more dissimilar 
to GR (DGP versus our simple tweak of the growth index) . They also exploited the ISW effect on CMB anisotropy, 
which we have ignored. Furthermore, we have included a recent prescription by Hu & Sawicki Q, which forces the 
matter power spectrum in MG to agree with GR on non- linear scales (see (jB13|) and preceding text); failure to account 
for this limiting behavior causes an over-estimation of the difference between the weak lensing observablcs in GR and 
MG. Each of these differences should lead to a less optimistic result in our case. 

In the next section, we introduce a new approach for comparing the 4 DES probes and show that this approach is 
much more successful at distinguishing our toy MG model from the canonical GR model. 



IV. MULTI-DIMENSIONAL CONSISTENCY TEST 



A. Searching for Parameter Tension 

Consider a simple example with one free parameter A and two probes. To determine if the probes are consistent 
we form 



X^(A)=E(A-A''OT:i(A-A«) (9) 



i=l 

where A*^*^ is the best fit value of the parameter from an analysis of probe i, and cr*^'^ is the error on probe i. A simple 
way to see if these two probes are consistent is to minimize the with respect to A. The value of A at this minimum 
is then the best fit value, and the value of there quantifies the goodness of fit in the standard way. There are two 
probes and one parameter corresponding to one degree of freedom, v = 1. So a Xmin ~ ^^y^ means that the fit 
is quite bad because the expectation value is (Xmin) = ior 1/ degrees of freedom. The difference, Ax^ = 9, tells us 
that the probes are inconsistent with 99.7% confidence. 

How can we forecast the expected tension between the 2 probes if the assumed model is incorrect? In that case, it 
will not necessarily be true that (A*^^^) = (A*-^^). Then the usual result, (Xmm) — ^! ^ill hold: the expected value 
of Xmin ^il^ larger, corresponding to a worse fit: (Xmin) = v + B with S > 0. A signature then of an incorrect 
assumed model is a value of B large relative to v. I.e., the large Xmin nicans that the probes are inconsistent in 
their measurements of A. For example, if i? = 9 and = 1, we would conclude that the 2 probes are inconsistent 
with 99.7% confidence. Thus, if we wish to forecast the tension among several probes, the tension parameter B is an 
appropriate quantity to calculate. 

Let us now generalize this l-paramcter/2-probc example to a situation where there arc M parameters and 
probes. Specifically, we are considering 8 cosmological parameters and 5 probes - the four of DES plus one from 
Planck. Suppose that probe i returns a best fit set of parameters Aq^*' with a covariancc matrix C^*^, and for now 
we assume that each of these is invertible so none of the probes is plagued with any degeneracy in parameter space 
(we will relax this assumption later). Let Aq, be a random point in cosmological parameter space. Then, a simple 
statistic [l3| which speaks to the agreement of the probes is 

x'(Aa) = E E(Ac - Ao^'^) (A^ - A^«). (10) 
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If the likelihood from each individual probe is Gaussian in parameter space and if the assumed model is correct, 
then the measured value of this will be drawn from a distribution with mean and variance equal to (TV — 1)M. 
Excessively large values of would falsify the underlying model. 



B. Predicting Parameter Tension 



We are now in a position to project how powerful a given set of probes will be when it comes to falsifying a model. 
Namely, we can compute the expectation value of as defined in Eq. ((TU)) if the true model were MG and determine 
by how much it exceeds (N — 1)M. First, we minimize and compute the value of A at its minimum: 



E 

/3,7 



(0 



(11) 



a/3 



We want to insert this into Eq. (|T0|) and take the expectation value. To do this, wc will set ((C^*-') ^) = i^'*', the 
Fisher matrix for probe i, and 



(12) 



where Aq'' = ^Aa*"*'^ is the expected outcome of the iih experiment. This assumes that the errors on the various 
probes are uncorrelated with each other, which is not strictly true since, for example, we expect a stronger weak 
lensing signal when we see more clusters, but it is a reasonable approximation (l9l [20j. Using these leads to 



(xL„) = (iv- i)M + ^^y„«y^W(F«)-i -^y„y,(G-i)„^ . 



i a/3 



where 



Va 



Ya = E^" 



(13) 

(14) 
(15) 
(16) 



If all the probes are expected to return the same parameter (the assumed model is correct) , then the two last terms 
on the right in Eq. ([T51) cancel, and (Xmin) = (-^ ~ l)Af . It makes sense therefore to subtract off a fiducial parameter 
set: 



AX — X Xfi'l 



so that 



(17) 

(18) 
(19) 



Then, the expected value of at the minimum depends only on these differences. In particular. 



with 



(Xmin) = (^-1)^^ + ^ 



i a./3 



(20) 



(21) 
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Our desired result - the projection for the excess due to inconsistency in the probes - has been reduced to 
a calculation of B. The ingredients of this calculation, Ay'^'\ are directly related to the errors in our theoretical 
predictions. By ([8]) and (fT8|) . we have 

Ay.«=^(Cov-i)«^AP,« (22) 

SO that Ay'*^ = Ay = when wc fit the correct model (or even an incorrect model as long as it does not produce 
tension in the probes) to the data. 

Using the probability distribution function for f degrees of freedom, we can use B to quantify the implied 
goodness-of-fit: we compute the probability of finding a worse Xmin than its expected value P{Xmin ^ (Xmin) i with 
1/ = Af TV — AI and (Xmin) = + B. As an example, with N ~ 5 and M = 8, we have v = 32. If we have chosen an 
incorrect model so that B = 14.2, then ^'(Xmin > (Xmin) ! ^ 0.046, meaning we predict that the constraints from 
the 5 probes will be inconsistent at the 95.4% or "2ct" level. 



C. Parameter Degeneracies 



In practice. Fisher matrices can be singular due to an experiment's insensitivity to a parameter or combination of 
parameters (parameter degeneracies). For example, because SN and BAO cannot measure non- geometric parameters 
(rif,. Us, (Tg), the corresponding rows/columns of the SN and BAO Fisher matrices are exactly zero (and the covariance 
matrix contains some corresponding infinite values). Also, the CMB cannot jointly constrain {wq, Wa, ^de, ^k) since 
these parameters essentially affect only one observable - the position of the 1st CMB peak - although there is some 
extra sensitivity to these parameters via the late ISW effect. Such insensitivity yields a singular and hence non- 
invertible Fisher matrix, or if the matrix is effectively singular (ill-conditioned), then its inverse formally exists but 
will be numerically unstable. These degeneracies represent more than a computational issue: they correspond to 
error ellipsoids which are infinite in some directions in parameter space, and hence parameter constraints cannot be 
inconsistent in those directions. E.g. we clearly do not expect the SN constraint to be inconsistent with any of the 
other probes in the erg direction. Subsequently, it is not sensible to count erg as a degree of freedom of the SN error 
ellipsoid. 

To handle degenerate Fisher matrices, we "clean" them using singular value decomposition. For each Fisher matrix, 
we find a unitary matrix U such that 

F = U'^AU (23) 

where A is a diagonal matrix whose diagonal elements are the eigenvalues of F. We replace the smallest elements of A 
with zeros (some are already zero) and recompute F using the above equation; this procedure negligibly changes the 
elements of F provided we only remove eigenvalues which are much smaller than the largest eigenvalue. We compute 
the inverse of F via 

F-'^ = U'^A-^U (24) 

where A"-'^ is diagonal and we have set 

^.-iN _ / l/Aca forAaaT^O , . 

J"" - j for A„„ = • ^^^> 

Let be the number of eigenvalues of i^^*^ which are zero (the nullity of the matrix). Starting again from (|13p . one 
now finds that 

{xLn)^iN~l)M-J2S^-'^+B (26) 

i 

with B unchanged from (|21l) . Effectively, the number of degrees of freedom v has been reduced by the total number of 
parameters that each probe cannot constrain. When very small eigenvalues of each F^^^ are set to zero, we find that 
B changes negligibly. This means that no significant tension is expected among the probes in these highly degenerate 
directions in parameter space, so we are indeed wrong to count them among the degrees of freedom. By setting the 
eigenvalues to zero, we force ourselves to evaluate tension only among the degrees of freedom where we expect tension. 
We excise the 3 most serious degeneracies from both the CMB and WL Fisher matrices, and we excise the 4 most 
serious degeneracies from the SN, BAO and CL Fisher matrices. Thus, when all of these probes are combined, the 
total number of degrees of freedom is = 5 x 8 — 8 — 18 = 14. 
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D. Results 

We suppose that we will mistakenly fit an S-parameter ACDM model to data in a Universe described by our toy 
modified gravity model. The MG model has a scale-independent linear growth history given by 7 = 0.68, and all 
ACDM parameters of the model are otherwise identical. In Table HI we show our estimates of the subsequent tension 
among the DES probes and Planck, assuming only statistical errors in the DES probes. As described in mV CI when 
a probe's individual Fisher matrix contains a significant parameter degeneracy, we do not count the corresponding 
direction in parameter space as a source of tension or a degree of freedom. 
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TABLE I: Expected tension between dark energy probes, assuming that our toy MG model is correct but we fit our GR+dark 
energy model to the data. We try different probe combinations - check-marks denote included probes. The total degrees of 
freedom, v, counts the number of non-degenerate parameter combinations measured by each probe, minus the one 8-parameter 
set used to minimize x^- The expected parameter tension, B, is calculated in (|2ip and vanishes when all probes are expected to 
have the same best-fit parameter set. We expect that (xmin) = + B, and this too-high value of Xmin ^ih cause us to interpret 
the constraints to be inconsistent (non-overlapping) when B is large. -P(Xmin > + B;u) is the expected goodness-of-fit of the 
probes, i.e. the probability that our probes, in a truly GR universe, would yield parameter constraints with more tension than 
the tension we predict due to fitting an incorrect model. 

The last row of Table H] shows that when all DES probes are combined with Planck, we expect their parameter 
constraints to have a goodness-of-fit of 0.0087. That is, the overlap of the error contours in the 8D parameter space 
will be so poor, we would interpret them to be inconsistent at about the 99.1% level. Using only WL, CL and CMB, 
we still expect these probes to disagree at about the 95% level. We predict no significant tension among the probes 
when either WL or CL is excluded. Our result is an improvement over the method of looking for overlap in the wq — Wa 
plane (see Figure[2]) since it leverages tension among all parameters, and docs not use the Planck prior multiple times. 
We remind the reader that our forecasts assume only statistical errors in the DES probes. Accounting for systematic 
errors would degrade parameter constraints, leading to more overlap, i.e. less tension. At the same time, we have 
chosen a modest toy MG model which differs from ACDM only in the linear growth of perturbations. Other MG 
models could easily produce more tension. 

Ours is an encouraging result in the search for modified gravity: it means that an inconsistency among the DES 
probes is, in principle, a useful diagnostic for identifying a wrong cosmological model. Of course, if such an incon- 
sistency were found, each probe's working group would revisit its pipeline, looking for systematic errors. If tension 
persisted, then the evidence for inconsistency would be strengthened. 

V. CONCLUSIONS 

The consistency of different dark energy probes promises to be a powerful tool in the quest to distinguish dark 
energy from modified gravity. Here we have illustrated that, using the Multi-dimensional Consistency Test (MCT), 
future probes from the Dark Energy Survey will be able to rule out standard (GR+dark energy) if the true gravity 
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model is only a modest modification of GR. Carrying out the MCT once the data are in reduces to computing the ^ 
of Eq. ([TU| while properly accounting for degeneracies as described in CI Although we have not explored this in 
detail here, projections of the MCT might make a useful metric for future surveys when trying to understand their 
constraining power towards modified gravity models, complementary to the figures of merit for dark energy [2l| . 
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Appendix A: Fisher Matrices for Unmodified Probes 
1. Supernova 

Given the number of supernovae N and a projected uncertainty on the magnitude of each supernova, we can form 
the Fisher matrix, 



N 

E 



d^i dfij 1 



(Al) 



where is the distance modulus for a supernova at redshift Zi and Pa are the parameters of interest. If we bin the 
supernovae in redshift bins of width Az = 0.1, each with an expected Na such that 



(A2) 



then the Fisher matrix becomes: 



dXa dXp al 



(A3) 



where the first bin is at z = 0.05 and the last at z = 1.25. We set a = 0.12 in low redshift (z < 0.8) bins and a ~ 0.24 
in the higher bins. Following the DES team, we also introduce a systematic floor, so that the error in any bin is 

t^bin = f^ll^a + Crioor ^ith (Tfloor = 0.02. 

The derivative of the distance modulus to bin a is 

d^ig _ 5 d\\idL{zg) 
dK ^ In(lO) dTc ■ 

For small values of the curvature, the luminosity distance can be approximated as 



dL{z) ~ (1 + z) 



6 



where 



with 



X{z) = 3000 Mpc 



dz 



(A4) 



(A5) 



(A6) 



E{z) EE h 



+ wq + Wg- — -) > + rtk{i + zY 

Ijo1 + ^ 1 + zJ 



1/2 



+ zf + Vtx exp <^ 3(1 + W0 + Wa) ln(l + z) - 3waz/(l + z)\ + 17fc(l + z) 



1/2 



(A7) 



We do not observe distance moduli directly but rather the peak apparent magnitude of the supernovae in a particular 
band, m{z) = Mq + ^(z), where Mq is a supernova's absolute magnitude. To account for uncertainty in the absolute 
magnitudes of type lA supernovae, we include Mq as a nuisance parameter in the SN Fisher matrix; Mq is completely 
degenerate with the Hubble constant h. The SN matrix can be combined with other matrices by padding the others 
with a row and column of zeros, adding all matrices together, and then marginalizing over Mg. 
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2. Baryon Acoustic Oscillations 

Galaxy surveys can measure a power spectrum of the galaxy distribution P{k) in a thin redshift bin around k — k^ 
with an expected covariance matrix 



where 



_^ P(fc.„)P(fcn) , 
Vn Veff 



(A8) 



Vm 



n{v)P{k) 
1 + n(r)P(fc) 



nP{k)R{n) 



nP{k)R{n) + 1 



Vo 



(A9) 



is an effective volume utilized for the measurement, Vq is the total volume of the survey, R = {1 + 13 fj?')'^ describes 
the effect of linear redshift space distortions, and 14 = is the volume of a shell in the Fourier space [for details 
see, [2^. Equation (|A9p holds if the number density of galaxies is constant within the volume. 



If we have a model of how the power spectrum depends on the cosmological parameters of interest Ai and if we 
assume that the errors on the power spectrum measurements are uncorrelated random Gaussian variables, then by 
propagating correlated errors given by the covariance matrix in Eq. (jASp we can get a Fisher matrix of cosmological 
parameters 



d^k dP{k)dP{k) 
(2^)2 d\, dX~ 



(AlO) 



[for details see [l^l • We want to forecast cosmological constraints from a measured position of the baryon acoustic 
peak in the three dimensional power spectrum. We will follow the formalism of Seo and Eisenstein [2^ [2^ . 
The baryonic part of the power spectrum can be modeled as 



P.(^)-^exp 



k 



(All) 



where sq is the sound horizon at recombination and fcgiik is the Silk damping scale. The Silk damping scale can be 
accurately fit by 



i.6(f}b/i')°''(an)°''[i + (io.4a„/i')-°°'']/i-' (/iMpc-i) 



(A12) 



As in |2a | we will multiply Eq. (jAlip by additional Gaussian functions to account for the erasure of information due 
to nonlinear evolution and photometric redshifts. The final Pb is given by 



sm 



Pb(fc) 



X cxp 



11*11 



1/2- 



1/2 



exp 



«silk 



1.4 



(A13) 



where s\\ and si_ are sound horizon scales measured along and across the line of sight, ji is cosine of the along the line 
of sight, P0.2 is galaxy power spectrum at = 0.2/iMpc~^, Aq is a normalization factor, and T,±_ model the loss 
of information due to nonlinear growth and Ez models the loss of information in radial direction due to photometric 
redshifts. As in [2^ we will use numerical values S|| = EoC?(l + /) and = SqG, where G is the growth function, 
/ = dlnG/dlna and Eg = 11.0/i~^Mpc for the cosmology with erg = 0.8. When the physical value of the sound 
horizon is known to high precision from CMB measurements the errors on s\\ and s± are equivalent to the errors on 
angular and radial distances Da and H . 

Derivatives of Eq. (jA13p with respect to sm and s± are 
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1^ ^ ^t^/„W, (A14) 

alns|| o\nx ' 

5Pb(a;) _ 9Pb(.T 



9 In s 9 In ; 

1/2 



(A15) 



where x = Ik'^^s'^^ + k'^s']_\ , /|| = /i^ and ,f±{lJ.) = 1 — /-t^. For most wavenumbers of interest kso is large and the 

sinusoidal terms oscillate rapidly. We will use this fact to replace cos^(fcso) term by its rms value of 1/2 and to drop 
terms proportional to sin^(fcso) when computing derivatives in Eqs. (jAl4p - (jA15p . Using Eq. ()A10|) with derivatives 
in Eqs. (|A14[) ~ (|A15|) results in a Fisher matrix on radial and angular distances H{z) and Da{z) 







X / fc^dfc 



exp h2(fc/fc3iik)i-1 







X exp 



[Pik)/Po.2 + {nPo.2Rit^))-'] 



k^{l - - k^^?Y.l exp [-fcV^S^] (A16) 



The Fisher matrix for H{z) and Da{z) computed from E q. (IA16P can be easily transformed into a Fisher matrix on 
cosmological parameters ujq, lUa, i^m, ^^A and h [see, e.g.. [27j. 



CMB 



Throughout, we compute Fisher matrices for the parameter set {wq, Wa, J^dE: h, ilb, rig, a^}- However, our prior 
from Planck CMB measurements comes in the form {wo,Wa,^T:)E,^k,^mh^,^bh^,nsA^P} where InP is the log of 



the amplitude of primordial fluctuations |16j . We transform the prior matrix to our preferred format using a simple 
chain rule. Starting from ([7]), 

- l^l^iC (A18) 
FLl^^ (A19) 



> 



(A20) 

Thus we express one set of parameters, A^, in terms of the others, A^, and compute the necessary matrix of partial 
derivatives. 

Appendix B: Fisher Matrices for Modified Probes 
1. Clusters 

DES will optically detect thousands of galaxy clusters out to a redshift of 1.6. The clusters will also be detected in 
microwaves by the South Pole Telescope (SPT) via the Sunyaev-Zeldovicli (SZ) effect, allowing tighter constraints on 
the cluster masses. Binning the detected clusters by redshift, the observables for our cluster experiment will be the 
total number of clusters in each bin, above a given mass threshold which allows detection by SPT. 

To conform with DES predictions, we assume the cluster mass function of Jenkins et al. psf with the parameters 
in their equation (B4). Let n{M, z) be the comoving number density of clusters with mass Al at redshift z. We take 

/Af \ Pco din CM J,, . , X 

= -M ThTM^^^'^) 
/(M,z) = 0.316exp(-l|log(aM(2)"')+0.67|3-82) , (B2) 
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FIG. 3: The SPT+DES number of detectable clusters Ni in redshift bins of width Az = 0.1 along with projected error bars in 
each bin. Solid curve assumes GR 7 = 0.55, and the dashed curve assumes modified gravity with 7 — 0.68. We plot the ratio 
of the two curves in the lower panel. 



Here, pco is the critical density today and cta/ is the RMS of the matter density field, smoothed by a top-hat filter of 
radius R where = SM/AiTpco- The RMS is calculated using linear theory: 

aj,{zf = I d\nk :^!i^l^4Mj2(fci?) ^ (B3) 
J 27r^ 

where Ji is the Bessel function of the first kind and Piin{k; z) is the linear matter power spectrum computed from 
the fitting formula of Eisenstein and Hu [2^. The minimum mass limit of clusters detectable by SPT, AI\i„j{z), was 
calculated in (30| . Let Ni denote the total number of clusters above this mass limit in the ith. redshift bin. It is given 

by 

N,='iTrf^y I dz^^;^ dMn{M,z) (B4) 

'Mii„(2) 



H{z) 



where x i^ comoving distance, Zi denotes the lower edge of the ith bin, and f^^^ = 0.125 is the sky coverage of the 
overlapping DES+SPT survey. We show Ni as a function of redshift for 7 = 0.55 and 7 — 0.68 in Figure [3] The 
clusters are divided into 16 redshift bins of width Az = 0.1, assuming perfect measurements of their redshift and 
mass. 

Our toy MG model only differs from ACDM via the linear density growth, parametrized by 7 in ([5]). In this case, 
computing the effect of MG on clusters is easily implemented by computing the linear growth function as a function 
of z for 7 = 0.68 and then using it to normalize Piin(fc;z) in (jBSp . A more general MG model could change the 
dynamics of collapsing halos, such as the halo formation time or the critical ovcrdensity for halo collapse [3l|. Studies 



of DGP [32|, [33| and /(i?)-gravity [3j| have shown that these models alter the critical spherical ovcrdensity, 5c^ by 
only 1-2% relative to ACDM. Furthermore, changes in halo formation times are already incorporated into the GR 
spherical collapse mass function of Sheth and Tormen [3l[, which has been shown to fit simulations well [e.g. Issl - iSTj . 
Therefore our assumption, that cluster numbers depend primarily on the linear growth function, is realistic. 

Assuming that the error on the number of clusters in the ith redshift bin is dominated by counting error, the 
covariance between bins is 

Cov[iV„ TVj] = %iV, . (B5) 
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The resulting projections in the dark energy parameter plane are shown in Figs. 2 when the underlying model is 
ACDM. Fig. 5 shows the projected constraints if the underlying model is MG with 7 = 0.68. Fig. 5 shows that to 
get the extra clusters that this MG model would produce, the fitter favors a larger w, thereby introducing more dark 
energy at early times. This tends to leave CTg closer to its present value, so mimics the increase in cluster abundance 
at high z. 



2. Lensing 

We begin with a review of relevant weak lensing equations [for a full introduction see e.g. [38l - l4l| . The lensing 
convergence at a particular sky position, can be expressed as the matter density contrast, S{x), projected over 

comoving distance, x, along the line of sight: 

/•OO 

/ dx6iex,x)W,{x). (B6) 
Jo 

The subscript i denotes the redshift bin from which source galaxies have been selected, and Wi{x) is the lensing kernel 
for that bin, defined below. The cosmic convergence power spectra and cross spectra, Ci-ij, are defined as 

= {27t)H^{1 + I') Ci.,j (B7) 

where I is the Fourier conjugate to 6, and we are working in the small angle limit, rather than decomposing the 
fields into spherical harmonics. Note that 5^ is a 2-dimensional Dirac delta function and the angle brackets denote 
an ensemble average. 

Measurements of galaxy ellipticities allow us to estimate the cosmic shear power and cross spectra, not the con- 
vergence, but these are equal at first order in the gravitational potentials. For N redshift bins, the weak lensing 
observables are the N{N + l)/2 total spectra for a given /. The leading order calculation of Ci-jj is 



dx_ 

dAixy 



Ci^^j^i -TT:T^w^(x)wAx)Ps{k-.x) (B8) 



where the function dA modifies distances in a curved Universe. For K = {H'^Q.k) ^ , 

X for K = (flat) 

dA{x)={ |if|-i/2sin(|i^|i/2;^) for A- >0 (closed) . (B9) 
\K\~^/'^smh{\K\^/'^x) for /v < (open) 

Ps {k; x) is the 3D matter power spectrum for k = £/dA{x) a distance X: accounting for the growth of structure. 
Equation (|B8p uses the Limber approximation, which assumes that the only matter density modes S{k) contributing 
to the lensing signal are those modes with k transverse to the line of sight. 
The lensing kernel is given by 

with Wo = §rim-ffo- Here, Pi{z) is the true (spectroscopic) distribution of galaxies in the ith redshift bin, and n^'^^i 
is the total projected number density of galaxies in that bin. Binning is done according to the galaxies' photometric 
redshifts (or "photo-z"s), which we take to be unbiased estimators of the true redshifts with gaussian scatter cr^. 
To avoid degrading parameter constraints, it is actually more important to have a small uncertainty in the bias and 
scatter of the photo-zs rather than a small scatter [42j . 

The redshift distribution of source galaxies in the survey is taken to be 

^(xz2exph(z/2o)'-']. (Bll) 

The median redshift of this distribution is Zmcd ~ ^0^/2, and the distribution is normalized so that the total projected 
number density of galaxies is 

dz ^ = V nS-\, = nS^i . (B12) 
dz 
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For DES, we take z,„od = 0.68, nS'*' = 12, and Uz = 0.08. The galaxies are divided into 7 photo-z bins with bin 
boundaries defined by z=(0, 0.37, 0.51, 0.63, 0.75, 0.89, 1.10, 2.0), which give them approximately equal projected 
number densities. 

Equation ^ describes the relation between the lensing potential $ — ^E* and the matter density 5 in a modified 
gravity scenario. Since our modified gravity model has / = 0, this relation is unchanged from GR, and the weak 
lensing power spectrum can be computed from the matter power spectrum via (jBSp . We need to compute the non- 
linear matter power spectrum in both the GR and MG cases. To compute the shape of the linear power spectrum 
in GR, we use the fitting formula of Eisenstein and Hu (29| . The redshift dependence of the linear power is given by 
the growth function, which can be easily adapted to incorporate the growth index, 7, in the MG case. We compute 
the nonlinear matter power spectrum in GR via the Halof it fitting formula of Smith et al. [43| ^, which takes the 
linear power spectrum as its input. Solar system observations indicate that that GR is valid on small non-linear 
scales, therefore a viable MG theory must agree with GR in that limit [13, and references therein]. Halof it does not 
impose such behavior on the non-linear power spectrum, therefore it must be adapted in the MG case. One way to 
accomplish this is to compute the linear power spectra in the GR and MG cases and then apply Halof it to both; 
a realistic non-linear spectrum for MG can then be constructed by interpolating between the two results so that we 
enforce the GR prediction for density modes that have gone very non-linear. Such an interpolation was proposed by 

i: 

^) = (B13) 

where Pgr is the non-linear power in GR and Pmgx comes from the naive application of a non-linear fitting formula to 
the linear MG power spectrum. The interpolation is moderated by two functions obtained from perturbation theory 



fc3 



i:\k,z)=^—Py^{k,z)^ , c„i(z)=A(l + z)"^ (B14) 

Here, P]\n(k, z) is the linear power spectrum in MG and the parameters ai, a2, and A are calibrated by N-body 
simulations (33, H^. For our working model, we take A = 0.3, ai = 1 and a2 = 0.16, which are their values in 
DGP; ignoring their weak dependence on 17 and erg. We remind the reader that we are not using DGP specifically, 
particularly since we are keeping the background expansion of ACDM. The ability of future weak lensing surveys to 
measure these nonlinear parameters was investigated by \at\ . 

The differences between the GR and MG predictions for Ci-^a are plotted in Figure |4] for a few redshift bins. These 
differences are the Ai-^ defined above, although for clarity, here we have divided them by the GR prediction. As 
expected, there is more power in the MG gravity model at early times. We impose a conservative cutoff of / < 1000 
since smaller scales contain non-linear and baryonic effects which are not completely understood even in GR, and we 
don't want incorrect predictions on these scales to be misinterpreted as violations of GR |4^ - l5H . 

To compute the weak lensing Fisher matrix, we need to compute the covariances between the observable spectra. 
The total observed power in a given redshift bin will be a combination of signal and noise: 

= (B15) 
where 7rms is the intrinsic scatter of one polarization of the galaxy shears. The covariance between the observables is 

r/^obs yoobs l //^obs i^obs i y^obsy^obs \ /"R1^^^ 

where AZ is the l-hm. width [5^ . Our results are insensitive to since the cosmic shear power spectra are relatively 
featureless. For DES, P'^^ = 0.12 and 7i.nis = 0.16. Figure [2] shows the DES projections for lensing. If the true model 
were MG (right panel) , then the lensing constraint would shift to more early dark energy (higher w) to accommodate 
the slower growth of structure. 



^ The formula was only calibrated for flat ACDM models and open models. It has been extrapolated to closed models and assumed to 
approximately hold for to ^ — 1. 
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FIG. 4: Fractional difTerence between weak lensing power spectra predictions for GR and our MG model with 7 = 0.68. Since 
we fix as, the models diverge with increasing redshift. Agreement is enforced towards large I by the interpolation in (|B13p . 



